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For many new accelerator applications, superconducting radio frequency (SRF) systems are the 
enabling technology. In particular for CW applications, much effort is being expended to minimize 
the power dissipation (surface resistance) of niobium cavities. Starting in 2009, we suggested a 
means of reducing the residual resistance by performing a thermal cycle [ll], a procedure of warming 
m a cavity after initial cooldown to about 20 K and cooling it down again. In subsequent studies 
[31, this technique was used to manipulate the residual resistance by more than a factor of 2. It 
was postulated that thermocurrents during cooldown generate additional trapped magnetic flux 
that impacts the cavity quality factor. Here, we present a more extensive study that includes 
measurements of two additional passband modes and that confirms the effect. In this paper, we 
also discuss simulations that support the claim. While the layout of the cavity LHe tank system 
is cylindrically symmetric, we show that the temperature dependence of the material parameters 
results in a non-symmetric current distribution. Hence a significant amount of magnetic flux can be 
generated at the RF surface. 


I. INTRODUCTION 

Much SRF research and development is presently con¬ 
centrating on finding means to reduce the losses in the 
cavity for next-generation accelerator applications. It 
was found that the surface resistance Rs can be impacted 
by a number of different treatments. Rs consists of a 
temperature dependent BCS term i?BCS and a residual 
resistance i?res, both of which may be impacted by the 
treatment. For example, it was discovered that the an¬ 
nealing of high RRR niobium with titanium or nitrogen 
can significantly change Rs and e^ecially the dependence 
of Rs on the RF magnetic field [ 3 ,0 • 

We expand here on studies Qiina to minimize the 
additional residual term in the cavity systems introduced 
by thermal currents based on thermal cycling. During a 
thermal cycle the superconducting cavity is warmed up to 
temperatures around 20 K which exceeds the transition 
temperature T^. Afterwards the cavity is cooled down 
again to the operating temperature. Depending on the 
temperature difference between the two ends of the cav¬ 
ity tank system during the superconducting (sc) phase 
transition, significant variations of the surface resistance 
are observed. In the best case the residual resistance of 
the 1/9 TT mode is decreased from 10.6 nD to 1.4 nil which 
is more than a factor of seven. 

In contrast to annealing techniques, the cycling proce¬ 
dure should not change Rbcs because the material prop¬ 
erties such as mean free path and coherence length are 
not influenced. Hence we attribute the observed modifi¬ 
cation of Rs for different cycles to a change in residual 
resistance. This assumption is borne out by extracting 
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the residual resistance from quality factor versus temper- 
ature measurements (see Section HTCl) . 


A. Thermoelectric effect in the cavity tank system 

The TESLA-type cavity reported on here received a 
heavy BCP (about 150 qm) prior to a 2 h bakeout at 
800°C. A light BCP etch followed the heat treatment. 
Before the helium tank was welded onto the cavity a 
quality factor of about 2 • 10^° in the tt mode at 2K 
was measured in a vertical test which corresponds to a 
residual resistance of « 1.2 nH if one assumes that Rbcs 
did not change between vertical and horizontal test (fit¬ 
ting parameters for Rbcs in horizontal test are listed in 
Section lll Cl) . In all tests presented in this paper, the cav¬ 
ity was installed horizontally in the HoBiCaT facility Q , 
similar to an installation in a cryomodule. 

In a previous study 0 , we suggested that thermoelec- 
trically induced currents and their associated magnetic 
flux is responsible for the change of Rs upon thermal 
cycling. It is dependent on the temperature difference 
along the cavity helium tank system during cooldown. 
In the horizontal setup, the system is fabricated of two 
materials: Niobium (cavity) and titanium (helium tank). 
If a temperature gradient is applied along the system, a 
voltage arises which drives a current along the cavity and 
back through the tank: 

U = (5Nb - ^Ti) ■ (Ti - T 2 ) (1) 

— ^ ‘ .^system, DC (2) 

where S is the thermopower or Seebeck coefficient of the 
respective material, Ti ^2 are the temperatures of the two 
joints, I is the thermocurrent and Rsystem, DC is the DC 
resistance of the system which will be dominated by the 
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resistivity of the helium tank, especially once part of the 
cavity goes superconducting. 

The magnetic flux associated with the current can then 
be trapped in the niobium when it becomes supercon¬ 
ducting and hence produces additional RF losses iSl. 

A number of questions were left open in our previous 
paper. For one, the thermopower for niobium and ti¬ 
tanium were poorly known so reliable estimates of the 
thermocurrents near were not possible. Furthermore, 
it was not clear how magnetic flux at the RF surface can 
be generated by a cavity tank system whose geometry 
is cylindrically symmetric. Finally, a systematic study 
of the surface resistance versus (controlled) temperature 
difference during the superconducting transition was still 
lacking. 

Here we present a far more systematic and in-depth 
study that addresses these questions. In the first part 
(Section |lTl) we discuss more extensive measurements of 
the quality factor of a TESLA cavity for many tempera¬ 
ture differences. Heaters were installed to generate tem¬ 
perature differences between the cavity ends that are typ¬ 
ically encountered during cooldowns in HoBiCaT. These 
measurements indeed confirmed a direct correlation be¬ 
tween the residual resistance and the temperature dif¬ 
ference (see Section umi). We also demonstrated that a 
reduced temperature difference (10 — 20 K) between the 
cavity ends during the very first cooldown of the cavity to 
the sc state can produce high quality factors without an 
additional thermal cycle being required. These results 
were confirmed for several modes of the passband and 
suggest that Q factors up to 10^^ are possible in our test 
cryostat, especially for the center cells where the external 
magnetic shield is most effective. 

In the second part (Section IHII) we discuss a quantita¬ 
tive analysis of the system to determine if indeed the ther¬ 
moelectric effect can explain the observed cavity behav¬ 
ior. Here, we turned to thermoelectric simulations of the 
cavity with helium tank. For these, accurate values are 
required for the temperature dependent thermoelectric 
properties of both niobium and titanium. Measurements 
presented in Section IHIBI now provide more reliable val¬ 
ues of the thermopower that were used in simulations. 

The simulations allow one to calculate the thermo¬ 
electric currents in the cavity tank system. Mechani¬ 
cally, the system exhibits toroidal symmetry (ignoring 
misalignment issues). Assuming perfect rotational sym¬ 
metry one would then expect a magnetic field between 
cavity and tank in azimuthal direction. Its magnitude 
is zero at the inner cavity RF surface, rising inside the 
wall, having a maximum at the outer surface of the cav¬ 
ity and falling back to zero on the outer surface of the 
tank. Hence it should not influence the RF properties. 
However, what counts is the symmetry of the current 
flowing in the system which in fact can be broken due to 
the temperature dependence of the material properties, 
as discussed in Section IlH Al Due to the broken symme¬ 
try, magnetic flux actually exists at the RF surface and 
within the cavity volume that potentially can be trapped. 


Our simulations, whose results are presented in Section 
IIII Cl confirm that the expected thermocurrent distribu¬ 
tion is consistent with the observed variations in cavity 
quality factor. 


B. Trapping of ambient magnetic field 

While we focus on the impact of thermocurrents and 
explain why it is a plausible explanation for the observed 
variation of Rg in our experiments, there is an additional 
effect that was intensively investigated during the last 
months: The change of the efficacy of the Meissner effect 
depending on the cooling rate and temperature distribu¬ 
tion for slow, homogeneous cooling 0,113, EMI. This 
effect can therefore also play a role in determining the 
residual resistance in cavities. 

Since our instrumentation is limited to sensors outside 
the helium tank we cannot make any statement on the 
exact temperature distribution of the cavity. However, 
we believe that a change in flux trapping efficiency of 
the residual field has only a minor impact on our par¬ 
ticular results. Our cavity was shielded from the earth’s 
magnetic field by a double (warm/cold) shield combi¬ 
nation, and the ambient field at the cavity position, as 
measured at room temperature, is low (typically a fair 
amount less than 0.5 pT at the center cells, about 1 pT 
maximum in the end cells 0] ). The inner cold shield was 
manufactured from Cryoperm, so that the shielding effi¬ 
cacy should theoretically be even better during cold op¬ 
eration. We measured the permeability of the material 
[Tsil and actually found little temperature dependence, so 
that the warm results should in fact be a reasonable esti¬ 
mate of the residual field during cavity operation. This is 
borne out by the best residual surface resistances of only 
1.4 nH observed in the center cells of the cavity when the 
temperature difference during cooldown is small and the 
thermoelectric effect should play no role. 

Furthermore we stress that the cavity we tested was 
not nitrogen doped. While doped as well as undoped 
cavities trap a higher percentage of ambient field when 
the cooling rate is decreased below a certain value, the 
same amount of trapped flux causes considerably less RF 
losses in an undoped cavity M- Hence the difference be¬ 
tween trapping all of the ambient flux in the center cells 
(at most 0.5 pT) and not trapping any of the flux would 
lead to a decrease of at most 1.75 nH in Rg for the 1/9 tt 
mode using the approximation of 3.5 nH/ pT [T3|. This 
is considerably less than the observed change of 9.2 nH 
between thermal cycles. For the tt mode the contribution 
may be larger and this is one of the reasons we measured 
several modes of the passband. We show later (e.g., Fig¬ 
ure HI) that indeed we observe no significant impact of a 
changing Meissner efficacy on our measurement results. 
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FIG. 1. Setup of TESLA cavity inside HoBiCaT cryostat with 
positions of Cernox temperature sensors 1-6 and heaters. 


II. THERMAL CYCLING EXPERIMENTS 
A. Experimental setup and cooling procedures 


A fully equipped TESLA-type cavity welded into a ti¬ 
tanium tank and with a TTF-3 input coupler installed 
was mounted horizontally inside the HoBiCaT cryostat. 
The cavity was equipped with Cernox sensors on the he¬ 
lium vessel head and beam pipes near the Nb-Ti joints, as 
well as helium inlet and outlet. Furthermore two heaters 
were attached, one on each beam pipe. The setup includ¬ 
ing the helium supply is sketched in Figure [TJ 

HoBiCaT can cool the cavity with different schemes. 
The cryoplant fills the helium tank via two valves: The 
filling valve (FV) at the bottom left and the Joule- 
Thomson valve (JTV) which fills the tank via the 2- 
phase-pipe from the top right. We used three different 
cooling schemes, the initial cooldown, the thermal cycle 
and the parked cooldown, which will be explained in the 
next sections. 

During the initial cooldown (whose temperature pro¬ 
file is shown in Figure [2^), mainly the FV at the bottom 
of the cavity was used for filling the helium tank with 
4.2 K helium. The coldest spot was always the one closest 
to the FV. This lead to high spatial temperature gradi¬ 
ents along the cavity and from bottom to top. The JTV 
was 25 % open as well to cool down the heat exchanger 
but no liquid helium reached the cavity. 

For a subsequent thermal cycle (as in Figure [2Jd) 
the FV was closed and remained so during the whole 
process. As the first step the JT valve was closed and 
the heaters on both ends were used to evaporate the he¬ 
lium in the tank and create a temperature difference be¬ 
tween the cavity ends if desired. The targeted difference 
could be adjusted by varying the heater power. Values 
chosen were typical of those encountered during normal 
cooldowns. Cernox sensor 5 (purple) measured a jump 
in temperature once the tank was empty and the cavity 
was normal conducting. In addition, we monitored the 





FIG. 2. Temperatures measured by Cx 1 - Cx 5 during ini¬ 
tial cooldown (a), a thermal cycle (b) and parked cooldown 
(c). The dotted line indicates the transition temperature of 
9.2 K. The dashed line indicates the time when the first sen¬ 
sor dropped below transition temperature and a 5 min interval 
before that transition. 


temperature of the evaporated gas with Cernox sensor 6 
(black) which exhibited a characteristic rise in tempera¬ 
ture (sensor is not displayed in Figure [5] for clearness). 
Now, the JT valve was opened again to start gas flow 
through the system and to restart cooling while the FV 
remained closed. 

To pinpoint the time of the phase transition we mea¬ 
sured the bandwidth of the tt mode during one complete 
warm-up/cooldown cycle. This allowed us to correlate 
the sc transition with the temperatures shown by Cer¬ 
nox sensors. The first sensor at the ends read close to 
9.2 K within less than one minute of the time when the 
first cells of the cavity went superconducting. 

Note that cycling scheme may lead to a reversed spatial 
temperature difference compared to the initial cooldown 
if not balanced by the heaters because the cold gas and 
liquid started to pour from the top right. As soon as the 
first liquid entered the tank a drop in temperature was 


































































4 



FIG. 3. Field distributions of three passband modes of a 
TESLA cavity: tt mode, 1/9 tt mode and 8/9 tt mode com¬ 
puted with CST and measured with a bead pull setup using 
an A1 sphere. 


registered by Cernox sensor 5 as can be seen in Figure [5] 

During one test series we performed a “parked 
cooldown” from room temperature which combines 
properties of both the initial cooldown and the thermal 
cycle (Figure [2t). The cooling procedure of the initial 
cooldown was adapted to stop well before the sc phase 
transition. The cryoplant was balanced to maintain a 
constant temperature at the FV for 48 h. The set point 
was first set to 30 K and then continuously lowered to 
14 K during this period. Both valves, the FV and the 
JTV, were used. After all temperature sensors were 
clearly in equilibrium (no temperature change with time) 
the set point was further lowered towards 1.8 K and the 
cavity tank system transitioned with a small AT < 10K 
into the sc state. 

Altogether, we performed one initial cooldown from 
room temperature, followed by 11 cycles, one parked 
cooldown from room temperature and another two cy¬ 
cles. After each cycle we performed a measurement of 
quality factor Q versus temperature at a (tt mode) gra¬ 
dient of 4 MV/m to extract the residual resistance for all 
three modes. 

The coupling during the RF cold tests was in most 
cases close to critical and always between 0.3 and 3.5 
which results in a low error margin below 10 % [a. A 
three stub tuner was installed in front of the input an¬ 
tenna and used to ensure optimal coupling for each mode. 


B. Cavity modes 


For the evaluation of the residual resistance in the 
various cells of the cavity we measured the quality fac¬ 
tor in three different passband modes: The tt mode 
(1299MHz), the 8/9tt mode (1298MHz) and the l/97r 
mode (1274 MHz). Each measurement is therefore an av¬ 
erage over the cells which are exposed to the RF field in 
the respective mode {Rs = G/Qo)- The geometry factor 


is taken from CST simulations: 

= 271.2 H 
G8^/9 = 271.5 H 
G^/9 = 268.3 H. 

The TT mode exhibits a homogeneous field distribution 
along the cavity. The 1/9 tt mode has the maximum field 
in the center cells and low field in both end cells whereas 
the 8/9 TT mode has minimum field in the mid cell and 
maximum field in the end cells. 

Figure [3] shows a calculation of the field distribution in 
the cavity computed with CST microwave studio. The 
plot is normalized to the same stored energy. We vali¬ 
dated the computed field distribution with a bead pull 
measurement after all RF measurements were completed 
(Figure 131). 

The bead pull measurement confirmed the overall field 
distribution of the three passband modes. Hence, they 
were used to evaluate contributions to the surface resis¬ 
tance depending on the location. We were able to dis¬ 
tinguish between end cell region (dominant in 8/9 7r) and 
mid cell region (dominant in 1/9 tt) but we could not 
determine the cell number because all modes are mirror 
symmetric with respect to the cavity center plane in the 
axial direction. 

The RF measurements confirmed that i?BCS did not 
vary between the three modes (see Section Hi Cl) . In con¬ 
trast, the measured i?res was influenced by localized con¬ 
tributions like the spatial variation of HoBiCaTs ambient 
magnetic flux due to inhomogeneous shielding at the end 
cells and maybe material defects. A higher surface resis¬ 
tance in the end cells lead to a lower Q in the 8/9 tt and 
TT mode than in the 1/9 tt mode. 


C. Extraction of residual resistance 

The temperature dependent BCS contribution to the 
surface resistance Rbcs was separated from the ternper- 
ature independent residual contribution R^es using [l8|: 

Rs{T) = A ■ exp + Rres (3) 

with B = (4) 

Kb 

For all measured modes, the fit parameters A and B over¬ 
lapped within the error margin confirming that i?BCS re¬ 
mained constant for different cycles. Hence we deter¬ 
mined overall BCS-parameters of A = (31.8 ± 2.2) 
and B = (15.7 ± 0.2) K corresponding to U = (1.70 ± 
0.02) ksTc = (1.35 ± 0.02) meV. With these parameters 
we extracted the residual resistance from all measure¬ 
ments. Figure |4] shows two examples for BCS fits and 
i?res extraction: Cycle 8 (highest temperature difference 
during phase transition) and the parked cooldown (low¬ 
est temperature difference during phase transition). The 
error bars of the RF measurement are within the symbol 
size. 
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FIG. 4. Arrhenius plots of cycle 8 (highest AT) and the 
parked cooldown (lowest AT). The thick lines indicate the 
BCS fits with A = 31.8 pH and B = 15.7 K for the stated 
i?res. The thin lines indicate the range of fits that are covered 
by the uncertainty of A and B. Rs was calculated using Rs = 
G/Qo and = 271.211, Gg^/g = 271.511, G^/g = 268.311. 


FIG. 5. Residual resistance measured with the three passband 
modes as a function of the temperature difference along the 
cavity at the onset of sc phase transition. The lowest residual 
resistance was achieved in the center cells following the cycle 
1 and the parked cooldown (l.dnll^Qo > 10^^ at 1.5K). 


Abcs ~ 0.9 nH 


^res 


at 1.5 K 

TT mode 

8/9 TT mode 

1/9 TT mode 

Initial cooldown 

13.6 

14.2 

10.6 

Thermal cycle 8 
(high AT) 

13.3 

18.6 

8.0 

Parked cooldown 

6.4 

9.9 

1.4 

Thermal cycle 7 
(low AT) 

6.5 

9.9 

1.9 


TABLE I. Examples of measured residual resistances of the 
three passband modes. 


D. Correlation of Rres with cooling parameters 

Table U provides examples of the measured residual 
resistances. We observed that cooldowns with a low 
temperature difference along the cavity result in lowest 
residual resistances. For the 1/9 tt mode (field dominant 
in center cells) this amounted to a residual resistance 
of 1.4 nfl. Presumably the quality factor then is domi¬ 
nated by the ambient field in HoBiCaT which is about 
0.2 —0.5 p.T along the center cells causing a residual resis¬ 
tance of 0.7 — 1.75nB [I^. Other, non-flux contributions 
may also add to the residual resistance. 

Figure [5] displays the extracted residual resistances for 
all cooldowns and cycles as a function of temperature 
difference AT. This is defined as: 

A rri ^Cxl + 70x2 7 cx3 + 7Cx4 


The temperature difference between the temperature sen¬ 
sors on the two ends of the tank is relevant for the ther¬ 
mocontact voltages. This parameter is a measure for the 
induced thermoelectric voltage and calculated at the in¬ 
stance when the first of the four sensors drops below the 
transition temperature. 

Figure [5] shows that the residual resistance of the cav¬ 
ity decreases with AT. The curves for the three passband 
modes run parallel indicating that all nine cells are sim¬ 
ilarly affected. Thus the change in i?res is a global effect 
consistent with the thermocurrent model. 

Furthermore, we see that the 1/9 tt mode has a signif¬ 
icantly reduced surface resistance compared to the two 
other modes while the residual resistance of the 8/9 tt 
mode is elevated. 

We believe this is due to an increased ambient field in 
HoBiCaT near the end cells. Measurements of the mag¬ 
netic shielding in HoBiCaT yielded that the magnetic 
field inside the end parts of the shield is increased due to 
cut outs for coupler and geometry effects (up to about 

ipT) a. 

Note that we were limited to only four temperature 
sensors outside the tank which only provides an incom¬ 
plete picture of the temperature distribution. Two sen¬ 
sors were near heaters which results in an additional off¬ 
set in AT. For more reliable values more sensors along 
the setup are needed. 

Figure [6] compares the average surface resistance of the 
modes to that of the tt mode. A linear dependency is vis¬ 
ible for the cycles and the parked cooldown. The initial 
cooldown for the 8/9 tt does not fit into the linear curve 
which is already visible in Figure [5l This might be ex¬ 
plained by the fact that the cooling dynamics of initial 
cooldown are different from the other coolings. The effect 
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FIG. 6. Rres obtained for the S/Qtt and I/Qtt mode as a 
function of Rres obtained for the tt mode. Linear regression 
leaving out the initial cooldown results in: 

Rs/ 9 tt = 1.2 • + 1.8 nfl and 

Rl/9n ~ 1.0 ' Rn — 4.7 nfl. 


would be largest for the cells in the higher ambient mag¬ 
netic field region, consistent with the observation that 
the 8/9 TT mode differs significantly. 

Leaving out the data point of the initial cooldown, the 
two graphs in Figure [6] exhibit a linear behavior. The 
slopes are 1.2 (8/9 tt mode) and 1.0 (1/9 tt mode) which 
is close to 1 meaning that all cells are affected by the ther¬ 
mocurrent in a similar way as would be expected since 
the current is constant along the length of the cavity. 
The current and resulting trapped magnetic flux adds a 
constant residual resistance to each cell. 

The results in Figure [5] and Figure [S] also demonstrate 
that the efficacy of the Meissner effect in expelling the 
constant ambient field plays no major role in our ob¬ 
served changes in Rg. Otherwise the 8/9 tt mode should 
change its residual resistance more than the 1/9 tt mode. 
Then the three curves in Figure [5] would not run parallel 
and the slopes of the curves in Figure [S] would not be 
linear with a slope near 1. The exception is the 8/9 tt 
mode result after the initial cooldown (Ai?s ~ 5 nSl) 
which may be due to the fact that the cooling conditions 
during initial cooldown significantly differ from the condi¬ 
tions during the thermal cycles and the parked cooldown 
resulting in different Meissner efficacies. This would most 
affect the Q value of the 8/9 tt mode since the measure¬ 
ment averages over the high-ambient-field end cells. 

The intercepts in Figure [6] are -1-1.8 nfl (8/9 tt mode) 
and —4.7nfl (1/9 tt mode). They are a measure of the 
different residual resistance that is not caused by thermo- 
electrically induced magnetic flux but is always present. 
The values indicate that the end cells of the cavity have 
higher average surface resistance. Sources can be either 
the increased ambient magnetic field near the end cells 
in HoBiCaT as discussed above or particulate contami¬ 


FIG. 7. Simplified geometry for GOMSOL simulations. Di¬ 
mensions are in millimeters. 

nation (e.g. from the nearby input coupler). 

The 1/9 TT mode which has the least field in the end 
cells exhibits a residual resistance down to 1.4 nil in the 
best case. The tt mode with equal field distribution in 
each cell shows an intermediate i?res of 6.4 nil while the 
8/9 TT mode with maximum field in the end cells has the 
highest Rres of 9.9 nil. If the low level of magnetic field 
in the center of the shielding would extend to the end 
cells, the low surface resistance, corresponding to a qual¬ 
ity factor of more than 10^^ at 1.5 K, might be achieved 
in the tt mode as well. 

In conclusion we can say that we observed a significant 
contribution to the residual resistance that increases with 
the temperature difference between the ends of the cavity 
tank system that is installed in a horizontal test stand. 
Given a high-Q cavity with an overall surface resistance 
of order 10 nil at 2 K this effect could dominate the losses. 


III. SIMULATION OF THERMOELECTRIC 
INDUCED MAGNETIC FLUX 

A. Setup for simulations: Symmetry and 
symmetry breaking 

To investigate the experimental results further and to 
make a quantitative analysis of the thermocurrents we 
turned to numerical simulations of the system. In par¬ 
ticular it is important to understand the distribution of 
the currents and the magnetic flux. In our previous pa¬ 
per, we made a rough estimate of the induced current 
and the corresponding magnetic flux. Now, we want to 
use temperature dependent material properties to yield 
a more reliable number on the amount of flux which is 
to be expected on the RF surface. For that purpose we 
conducted COMSOL simulations which will we discuss 
in the next section. 
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FIG. 8. Different areas in which thermal boundaries could 
be applied in COMSOL simulations. No boundary condition 
was applied to the rest of the setup. 


We modeled a simplified geometry consisting of two 
coaxial cylinders (Figure [7]). The inner cylinder which 
resembled the cavity had the typical cavity wall thick¬ 
ness of 2.8 mm and was defined as niobium. The outer 
cylinder and the two end plates resembled the helium 
vessel and helium vessel heads and were defined as tita¬ 
nium with a thickness of 5 mm (vessel) or 10 mm (vessel 
heads) respectively. The 1 m length of both cylinders was 
the approximate length of a TESLA cavity in a tank. The 
model exhibited perfect cylindrical symmetry. 

The simulations were performed in several steps. For 
each step different temperature boundary conditions 
were applied. We defined ten areas on the outer surface 
of the tank (Figure|8]). Each area could either be assigned 
a hxed temperature for the simulation or left undefined 
as the rest of the setup. Thereby we created temperature 
distribution similar to our experimental situation. 

In the first step, we defined the temperatures only on 
the end plates, areas 1 (10K) and 6 (100 K). In this 
case the isothermals were perpendicular to the cylinder 
axis. The resulting current distribution was cylindrically 
symmetric and the azimuthal magnetic flux on the RF 
surface was zero though it increased rapidly as one moves 
into the wall. This situation is as expected from Gauss’ 
Law. 

In the next steps, the areas 2-5 and 7-10 were 
assigned fixed temperatures as well, leading to asymmet¬ 
ric temperature distributions. The isothermals were not 
orthogonal to the cylinder axis anymore but distorted. 
Thereby the electrical symmetry was broken due to the 
temperature dependence of the electrical resistance. This 
resembles the cavity cooldown when liquid helium started 
to collect first at the bottom of the tank. 

Simulations where the inner cylinder is replaced by the 
actual cavity geometry are not presented in this paper 
but are in preparation. 


B. Input data 

For these simulations the temperature dependent ma¬ 
terial properties are needed. We used literature data as 
given in references (Tol - f^ for the temperature dependent 
electrical conductivity. The thermal conductivity was set 
constant to improve the convergence of the COMSOL 
simulations. We choose kns = 50-^ and kti = 22-^ 
as approximations [20l - l22l| . 



FIG. 9. PPMS thermopower data of niobium and titanium 
for the small (2 mm x 2.3 mm x 30 mm) and the large sample 
(2 mm X 5 mm x 30 mm). 


On the thermopower only poor information is avail¬ 
able. Hence we performed sample measurements of the 
thermopower on niobium {RRR = 300) and titanium 
(grade 2). The dimension of the samples were approxi¬ 
mately 2mmx5mmx30mmand 2mmx2.3mmx30mm. 

The experiments were performed with a physical prop¬ 
erty measurement system by Quantum Design using 
the thermal transport option |24| . 

The system allowed high precision measurements of 
the thermopower over a wide temperature range. We 
measured all four samples (two niobium, two titanium) 
in the range from 4.5 K to 300 K several times. 

Figure [S] shows the measured thermopower S'[|j,V/K] 
as a function of tem per ature. The data is in agreement 
with literature data [251 - 1^ . For the titanium data, the 
small and the large sample were the same. The curve 
was fitted as shown in Figure [HI A table was extracted 
and entered into COMSOL. 

For niobium, the two samples did not give the same 
result. In general, both curves are similar and have the 
same zero-crossing but the amplitude was higher for the 
large sample compared to the small sample. The COM- 
























FIG. 10. Results of the COMSOL simulations. The left cylinders display the magnetic field distribution inside the cylinders at 
r = 102 mm. The legend is given on the left reaching from « 50nT (green) to « 13 p.T (red). The right figures display a cut 
through the temperature distribution inside the cylinders. The legend for the temperature reaches from 10 K to 100 K. 


SOL simulations were first performed with the fitted data 
for the small sample as shown in FigurelHlto obtain a con¬ 
servative lower limit. Note that the disagreement in the 
measured thermopower in the temperature range up to 
75 K is only in the 30 % range. The simulations were also 
repeated with the less conservative Seebeck values for a 
sensitivity check, but we found that the impact was not 
dramatic (see later Figure[T3]). Nevertheless we intend to 
repeat the measurements in the future to provide more 
accurate results. 


C. Magnetic flux at the RF surface 

We performed hve simulations: 

1. A symmetric test distribution as a consistency 
check for the simulations. Surface 1 was set to 10 K 
and surface 6 to 100 K. 

2. Asymmetric studies: Surfaces 1 and 10 at 10 K and 
surfaces 5 and 6 at 100 K 

3. Surfaces 1, 10 and 9 at 10 K and surfaces 4, 5 and 
6 at 100 K 

4. Surfaces 1, 10, 9 and 8 st 10 K and surfaces 3, 4, 5 
and 6 at 100 K 

5. Surfaces 1, 10, 9, 8 and 7 at 10 K and surfaces 2, 3, 
4, 5 and 6 at 100 K 

The simulations were performed as two steps. First, 
the temperature distribution and the resulting electric 


current distribution was computed. In the second study, 
this input data was used to calculate the magnetic field 
distribution. 

The RF surface is at a radius of r = 103.3 mm. Fig¬ 
ure [ini shows the solution for the magnetic field inside 
the cylinders close the RF surface at r = 102 mm to en¬ 
sure that it is inside the cylinder and not in the cylinder 
wall due to meshing. The arrows display the direction of 
the magnetic field at the same radius. The arrow size is 
proportional to the field strength. Generally the field is 
near perpendicular to the surface in high field regions. 

In addition to the magnetic field distribution. Figure 
[inidisplays cuts with the temperature distributions inside 
the cylinders. 

Figured!] displays the magnetic field at a cut through 
the walls of the two cylinders at z = 0.5 m. For the sym¬ 
metric case, no field is inside the inner cylinder, the field 
increases in the inner cylinder wall and decreases towards 
the outer cylinder. Outside of the two cylinders there is 
no magnetic field anymore. This is just as expected from 
Gauss’ Law. For studies 2 to 5, the increasing asymmetry 
of the current density leads to an increased magnetic field 
inside the cylinders. This field can get trapped during su¬ 
perconducting phase transition and produce RF losses as 
we have observed in the experiment. 

Figure [121 displays a cut of the RF surface along the 
z-axis. The position of the cut was chosen to be at 
X = 102 mm, y = 0 mm. Furthermore the Figure dis¬ 
plays the temperature difference between top and bottom 
of the inner cylinder. Where the temperature difference 
is maximum, the inhomogeneity in the electric current 
distribution is largest and hence the most magnetic field 
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FIG. 11. Absolute magnetic field along a cut trough the cylin¬ 
ders at a = 0.5 m. 



FIG. 12. Top Figure: Absolute magnetic field along « at 
X = 102mm and y — 0 mm. The peaks at z = 0 mm and 
z = Im are due to the sharp corners when the cylinder passes 
into the end plates. Bottom Figure: Temperature difference 
between top and bottom of the inner cylinder. 


penetrates the RF surface. 

The results clearly illustrate that a non symmetric tem¬ 
perature distribution leads to a nonzero magnetic field on 
the RF surface. The highest fields are obtained where the 
temperature difference between top and bottom part are 
largest. In study 5, this effect leads to a more or less 
constant magnetic flux along the cylinder while in study 
3 where the temperature gradient from bottom to top is 
less distinct the magnetic field increases towards the cen¬ 
ter of the cylinders {z = 0.5 m). In the experiment, the 
highest degree of broken symmetry is obtained when the 
first superconducting path through the cavity establishes 
below 9.2 K. This corresponds to an extreme example of 
case 5 where the DC resistance of the bottom part of the 
inner cylinder drops to zero while the rest of the cavity 
is still normal conducting and hence the thermocurrent 
is forced into the sc path. In this case all cells would be 
affected by the additional induced flux in the same way 
which is in agreement with the linear increase in Fig¬ 
ure El with slope 1. Temperature mapping data of a 
dressed multicell cavity in a horizontal test would give 
more insight in the real temperature distribution. 

Finally, we varied the temperature difference in the 
simulation and set the warmer boundaries of study 5 to 
the temperatures 80 K, 60 K, 40 K and 20 K while the 
other surfaces remained fixed at 10 K. The result is 
given in Figure [13] together with the experimental find¬ 
ings. Since the measurements of the thermopower of nio¬ 
bium allowed for 30% higher values, we added bars in¬ 
dicating the magnetic field at the RF surface for a 30 % 
higher thermopower of niobium. The results show that 
the magnetic field increases with the temperature differ¬ 
ence. While we do not claim that the quantitative values 
should be directly compared, the trend is very similar. 

In conclusion, we assess that the COMSOL simulations 
support the assumption that thermoelectrically induced 
magnetic field penetrates the RF surface due to an inho¬ 
mogeneous temperature distribution. We see that even 
for a moderate asymmetry the flux is locally of order 
several pT. The RF losses resulting from the increased 
magnetic field will depend on the fraction of the nine cell 
cavity that is affected by the thermocurrent due to more 
pronounced asymmetry, the exact location of the induced 
magnetic field relative to the RF field and the amount of 
flux that eventually gets trapped in the material. How¬ 
ever the variation of about 10 nD seen in the cavity tests 
are consistent with the thermocurrent hypothesis. 


IV. SUMMARY AND OUTLOOK 

We demonstrated that poorly controlled cooling con¬ 
ditions can significantly deteriorate the quality factor of 
SRF cavities due to the generation of thermocurrents. 
While already a simple misalignment of the cavity to the 
helium tank could break the symmetry of the cavity tank 
system, although the effect on Rs should be small [2^ 
even a mechanically cylindrically symmetric system will 
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FIG. 13. Magnetic field at the RF surface simulated for dif¬ 
ferent temperature differences in combination with the data 
of Figure [5] The red bars indicate the magnetic field that is 
obtained assuming a 30% higher thermopower, correspond¬ 
ing roughly to the uncertainty resulting from the discrepancy 
depicted in Figure [9] 


generate magnetic flux on the RF surface. The presented 
simulations explain how the electrical symmetry is bro¬ 
ken due to the temperature dependance of the electrical 
conductivity. 


We also demonstrated that an initial cooldown through 
Tc with small A T can in principle yield residual resis¬ 
tances at 1 nfl (parked cooldown), provided the external 
magnetic shielding is very effective in eliminating exter¬ 
nal residual flux. 

What we did not investigate in this paper is the influ¬ 
ence of changed cooling conditions on the efficacy of the 
Meissner effect as we explained in section im since our 
magnetic shield eliminated most of the residual field. 

Temperature mapping during dressed cavity horizon¬ 
tal test could improve the understanding on the real tem¬ 
perature distribution during cooldown and cycles. Fur¬ 
thermore it would help with characterization of cooling 
parameters like cooling speed, temperature gradients in 
all spatial directions and homogeneity. 
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